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Let Y be an Ornstein-Uhlenbeck diffusion governed by a sta¬ 
tionary and ergodic Markov jump process X: dYt = a(Xt)Yt dt + 
it (At) dWt, To = Vo- Ergodicity conditions for Y have been obtained. 

Here we investigate the tail propriety of the stationary distribution 
of this model. A characterization of either heavy or light tail case is 
established. The method is based on a renewal theorem for systems 
of equations with distributions on R. 


1. Introduction. The discrete-time models Y = ( Y n ,n £ N) governed by 
a switching process X = (X n ,n € N) fit well to the situations where an 
autonomous process X is responsible for the dynamic (or regime ) of Y. 
These models are parsimonious with regard to the number of parameters, 
and extend significantly the case of a single regime. Among them, the so- 
called Markov-switching ARMA models are popular in several application 
fields, for example, in econometric modeling [see Hamilton (1989, 1990)]. 
More recently, continuous-time versions of Markov-switching models have 
been proposed in Basak, Bisi and Ghosh (1996) and Guyon, Iovleff and Yao 
(2004), where ergodicity conditions are established. In this paper we investi¬ 
gate the tail property of the stationary distribution of this continuous-time 
process. One of the main results (Theorem 2) states that this model can 
provide heavy tails, which is one of the major features required in nonlinear 
time-series modeling. Note that heavy tails may also be obtained by using 
a Levy-driven Ornstein-Uhlenbeck (O.U.) process (without Markov switch¬ 
ing); see Barndorff-Nielsen and Shephard (2001) and Brockwell (2001). 

The considered process Y, called diffusion with Markov switching , is con¬ 
structed in two steps: 
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First, the switching process X = (X t )t> o is a Markov jump process [see 
Feller (1966)], defined on a probability space (fl,Al, Q), with a finite state 
space E = {1,..., N}, N > 1. We assume that the intensity function A of X is 
positive and the jump kernel q(i,j ) on E is irreducible and satisfies q[i, i) = 
0, for each i € E. The process X is ergodic and will be taken stationary with 
an invariant probability measure denoted by //. 

Second, let W = {Wt)t >o be a standard Brownian motion defined on a 
probability space (0, B , Q'), and let T = (Et) be the filtration of the motion. 
We will consider the product space (H x 0, A x B, ( Q x g Q')), P = Q ®Q' 
and E the associated expectation. Conditionally to X , Y = (Y t ) t >o is a real¬ 
valued diffusion process, defined, for each uj 6 17 by: 

1 . Yq is a random variable defined on (Q,B,Q'), JFo-measurable; 

2. Y is solution of the linear SDE 

(1) dY t = a(X t )Y t dt + a(X t )dW t , t> 0. 

Thus (If) is a linear diffusion driven by an “exogenous” jump process (Xt). 

We say a continuous- or discrete-time process S = (St)f>o is ergodic if 
there exists a probability measure m such that when t —> oo, the law of 
St converges weakly to m independently of the initial condition Sq. The 
distribution m is then the limit law of S. When S is a Markov process, m 
is its unique invariant law. 

In Guyon, Iovleff and Yao (2004), it is proved that the Markov-switching 
diffusion Y is ergodic under the condition 

(2) a = a{i)ii{i) < 0. 

i£E 

The main results of the present paper are the following theorems. Note that 
Condition 2 will be assumed satisfied throughout the paper and we denote 
by v the stationary (or limit) distribution of Y. 

Theorem 1 (Light tail case). If for all i, a(i ) < 0, then the stationary 
distribution v of the process Y has moments of all order; that is, for all s > 0 
we have 

/ \x\ s u(dx) < oo. 

JR 

Theorem 2 (Heavy tail case). If there is an i such that a(i ) > 0, one 
can find an exponent k > 0 and a constant L > 0 such that the stationary 
distribution v of the process Y satisfies 

f'V(]f,+oo[) —•> L, 

c—>+oo 


t K v(] — oo 


t —>+oo 
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Note that the two situations from Theorems 1 and 2 form a dichotomy. 
Moreover, the characteristic exponent k in the heavy tail case is completely 
determined as following. Let 


si = min 


A« 


a(i) 


a(i) > 0 


M s = (q(i,j) 


\{i) 


X(i) - sa(i)) ijeE 


for 0 < s < si. 


Then k is the unique s G]0, si[ such that the spectral radius of M s equals 
to 1. 

The proof of Theorem 1 is a consequence of a result of Guyon, Iovleff and Yao 
(2004), and the proof of Theorem 2 is based on a recent renewal theorem 
for systems of equations reported in de Saporta (2003) and on an AR(1) 
recurrence equation satisfied by the discretization of Y that we will define 
in Section 2. In Section 3 we study an operator related to our problem and 
prove Theorem 1. Sections 4-7 are devoted to the proof of Theorem 2. First 
we state two renewal theorems for systems of equations. Then in Section 5 we 
derive the renewal equations associated to our problem. In Sections 6 and 7 
we prove Theorem 2, the latter section being dedicated to the proof that the 
constant L is nonzero. Finally, in Section 8 we give further details on the 
computation of the exponent k. 


2. Discretization of the process and an AR(1) equation. First we give 
an explicit formula for the diffusion process. For 0 < s < t, let 

$(s,t) = $ Sj t(o;) =ex.pj a(X u )du. 

The process Y has the representation [see Karatzas and Shreve (1991)]: 


Y t = Y t (u) = <f>{0,t) 


Y 0 + f <5>{U,u)- l o{X u )dW u 
Jo 


and for 0 < s < t, Y satisfies the recursion equation 


Y t = ${s,t) 


Y s + 4>(s, u)~ 1 a(X u ) dW u 

= d*(s, t)Y s + f exp f a(X v )dv a(X u )dW u . 

J s L Ju 

It is useful to rewrite this recursion as 

(3) Y t (u) = 

where Gs.i is a standard Gaussian variable, function of (W u , s <u<t), and 


V s ,t( u ) = / ex P 2 / a(X v )dv 

J S . Ju 


a 2 (X u ) du. 
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For 5 > 0, we will call discretization at step size 6 of Y the discrete-time 
process Y^ = (Y n $) n , where n € N. Our study of Y is based on the investi¬ 
gations of these discretization Y^ as in Guyon, Iovleff and Yao (2004). 

More precisely, for a fixed 6 > 0, the discretization follows an AR(1) 
equation with random coefficients: 

(4) *(„ + l)«M = <bn+lMYn,5M + V^)i n+U 

with 


<f> n+ i(a;) = <f> n+ i(<5)(a;) = exp 


r(n+l)<5 


.JnS 


a(X u (u)) du 


r(n+l)S r r 

Vn+i (u) = exp 2 / 

JnS L Ju 


(n+l)8 

2 / a{X v {io))dv 


a 2 (X u (u)) du, 


where (£ n ) is a standard Gaussian i.i.d. sequence defined on {Q,B,Q'). Note 
that under Condition 2, all these discretizations are ergodic with the same 
limit distribution v [see Guyon, Iovleff and Yao (2004)]. 


3. Study of a related operator. We now introduce a related operator A 
and investigate its properties. Fix s > 0 and 5 > 0. We define the operator 
A ( s ,s) by 

A (s> 5 ) <^=E^(%(X 5 )], 

for every function ip : E —» R and every i in E. It has the following semigroup 
property: 


Proposition 1. Fix s > 0. Then for all 5 ,7 > 0 we have 

^■(s,5)^(s, 7) ^ ^-(s,<5+7) ■ 

Proof. Set 92 : E —>■ R and i in E. We have 

A (s,S) A (s, r y)F('i') 

= E i [$;(5)E Xi [$;(7)^ 7 )]] 

= Ej exp J a(X u ) dii'jExg exp J a(X u ) dvJj yp(Xy) 
Then the Markov property yields 

A {s,S) A {s,AT^S) 

= Ei[$f(<J + 7M^5 +7 )] 

A {s,8+ r y)T^S) • ^ 


rS +7 

exp [ s a(X u ) du )ip(Xs+^ 


Note that A^ s ^ip{i) = ®i[ < bflx 5 =j] ( ^(j); and therefore A( s ,S) can 

be rewritten as the matrix (( A ( s ,S))ij)l<i,j<N with {A^ S ))ij = Ej[$flx 5 =j]. 
Note also that it is a positive operator. 
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3.1. Spectral radius. Now we investigate the properties of the spectral 
radius of A. First, we recall a result from Guyon, Iovleff and Yao (2004). 

Proposition 2. Fix s > 0 and 5 > 0. Then Ar s> $-\ is irreducible, aperi¬ 
odic and satisfies 

(5) E^[($! • • • <t> k ) s ] = £ 4, 5 ) 1(*M0 = Mm)1, 

ieE 

where 1 is the constant function equal to 1 on E. 

We denote by p(X) the spectral radius of a matrix X. Proposition 2 yields 
the following corollaries. 


Corollary 1. We have 

P(A M )= lim (E li [(* 1 -* k ) a )) 1/k . 

k—> oo 


Proof. As A( Si , 5 ) is a (component-wise) positive, irreducible and aperi¬ 
odic matrix, Theorem 8.5.1 of Horn and Johnson (1985) gives the existence 
of a matrix Br s ^ with positive coefficients such that 


( 6 ) 


OW" : „ _ 

{p(A( StS ))) n n ^°° M) 


This result and (5) yield the expected result. □ 


Corollary 2. For all fixed 5 > 0, the mapping s i —> logp(.A( S( $)) is 
convex on M + . 


Note that for all fixed 5 > 0 and i in E, we have -A(o,<$)l(i) = Ej(l) = 1. 
Thus, as A(o,< 5 ) is a positive operator, it is also a stochastic matrix and 
P(^(0,<5)) = 1- 


Proposition 3. For all fixed 6 > 0, the right-hand derivative of the 
mapping s i—> logp(A( s at 0 is negative. 

Proof. As all the functions considered are convex, we have 

A r\ -t 

Tplog (p(A(s,S)))= lim w—logE^-Fi • • •<F n ) K ] 

OS n^oo Os n 

_ lim l E At [(4> 1 ---^.E'Lilog^] 

J^o n E M [(4>i---4> n ) K ] 
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The sequence (<3? n ) is stationary, thus the ergodic theorem yields 


(7) 


n 


E log—> E„[log$il, F u -almost surely. 

n.—>no ^ J ^ 


k=1 


But E ft [log<l>i] < 0 because of Condition 2. Indeed, we have 


E M [log3>i] =E m 


r8 1 rS 

/ a(X u )du = / Eu[a(X u )] du = 5a < 0. 
Jo . Jo 


Thus we get, as expected, 

d_ 

ds 


1 , 


log (p(^M))) = 11™ - e m 

s=0 n^oo n 


i°g ^ 


i=1 


= E At [log<hi] < 0. 


□ 


Corollary 3. Fix 6 > 0. We have the following dichotomy: 

(i) either for all s > 0, < 1, 

(ii) or there exists a unique k > 0 such that p(A^ fiS ' ) ) = 1, and in this case 
p(A( S ’, 5 )) > 1 for all s> k and p(A( s ^) < 1 for all 0 < s < k. 

3.2. Choice of 5. Now we are going to prove that the preceding di¬ 
chotomy is in fact independent of the value of 5. 


Proposition 4. Fix s > 0. The following propositions are equivalent: 

(i) there exists 5 > 0 such that p(Ar g ^) < l, 

(ii) for all 6 > 0 we have p(A/ s g^) < 1. 

The same equivalence is true if we replace “< 1” by “> 1” or “= 1.” 


Proof. Set 5 > 0 such that p ( A^ s ^)) < 1, and 7 > 0. For all integer 
n > 1 we define m n £ N* and 0 < /3 n < 6 by 717 = m n 5 + j3 n ( m n the integer 
part of 717/5 and (3 n its fractional part multiplied by 6). Thus Proposition 1 
yields 

7l(<s, 7 ) = A^ sn 7 ) = A( s ™ S j A( s fi n y 

But for all n we have 


II^VaoII <maxEi[$f(/? n )] 

< exp ^s/3 n max(aj) 

< exp I S(5max(aj) 
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This upper bound is independent of n. Thus we have 

log Pf S)7) II < log 11^)11+c, 
where c is a positive constant. We get 

logp(^( Si7) ) = lim-logP^ 7) || 

< lim sup ^ log 11^)1 

n ft v ’ ' 

= -logp(A( s>5) ), 

as m n ~ nyd” 1 . Hence p{A^ s ^) < p{A^ s ^) l/5 < 1. 

For the case “= 1,” fix So and a corresponding n such that p(A( K # g \) = 1. 
The mapping s \—> p(Au ^) is log-convex hence continuous. Thus we have 

P( A (k,So)) = su PP( A (s,8 0 ))- 

S<K 

Set 8 > 0. We want to prove that p(A^ K ^) = 1. According to Corollary 3, 
for all s < k we have p(A^ s ^ 0 )) < 1. Thus the preceding study yields that for 
all s < k we also have p(A( s ^) < 1. Hence we have 

p( A (n,S)) = sup p(A (Si<S )) < 1. 

S<K 

Suppose that p{A/ K ^) < 1; then the first case implies again that p(Au ^) < 
1, which is impossible. Thus we have p(A^ K = 1 as expected. 

The case “> 1” is a consequence of these two cases and Corollary 3. □ 

In the following we will write A s instead of each time it is nonam- 

biguous. We have an easy criterion to know in which case we are. 

Proposition 5. The following properties are equivalent: 

(i) for all i in E, a(i) < 0, 

(ii) for all s > 0 ; p(A s ) < 1. 

Proof. Suppose that for all i in E we have a(i) < 0. Fix 5 > 0. Then for 
all s > 0, we have <3?f < 1. Thus for all i, A s l(i) = Ej[<J>f] < 1, and component¬ 
wise we have A s 1 < 1, which implies that p(A s ) < 1 for all s > 0. Corollary 3 
then yields that for all s, we have actually p{A s ) < 1. 

Now suppose there exists an io such that a(*o) > 0. Fix s > 2A(*o)a(io)~ 1 - 
It is proved in Guyon, Iovleff and Yao (2004) that for all function p from E 
into M and all i in E we have for small 8, 

(8) A s p(i) = [1 + <5(sa(z) - A(i))]</?(z) + 8\{i) ^2[q(i,j)ip(j)\ + o(<5). 

3¥=i 
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Let if be the function from E into R such that if(i o) = 1 and if(i) = 0 for all 
i / io- Then for all i ^ io we have A s if(i) = Ej[$f lx 5 =i 0 ] > 0 and for i = io 
we have 

A s if>(i 0 ) = 1 + 5(sa(i 0 ) - A(*)) + o(5) > 1 + S sa ^ + o(<5) 

as we have chosen s > 2A(i)a(*o) _1 • Thus component-wise, for small enough 
5, we have 

Ag'ip > ^1 + 5 Sa ^ + o(<5)^ if) 


> 





Thus p(A s ) > 1 + > 1. □ 


This proposition ends the proof of Theorem 1 since we have the following 
result from Guyon, Iovleff and Yao (2004) that relates the spectral radius of 
A s to the moments of the stationary law v. 


Proposition 6 . Set s> 0. If p(A s ) < 1, then the stationary law u of Y 
has a moment of order s. 


The proof of Theorem 1 is now complete. 


4. Renewal theory for systems. Now we proceed to prove Theorem 2. 
From now on, we will assume that there is an i such that a(i) > 0. Our 
approach is based on a new renewal theorem for systems of renewal equa¬ 
tions. First we introduce some notation and conventions that we will apply 
throughout. 

Let F = be a matrix of distributions: nondecreasing, right- 

continuous functions on R into R+ with limit 0 at — oo. For all p x r matrix 
H of Borel-measurable, real-valued functions Hij on R that are bounded on 
compact intervals, we define the convolution product F * H by 


(F*H) ij (t) = J2 / H kj (t 


u)F ik (du) 


where it exists. 

The transpose of a vector or matrix X will always be denoted X'. We 
study the renewal equation Z = F * Z + G, where G = (G i,..., G p )' is a 
vector of Borel-measurable, real-valued functions, bounded on compact in¬ 
tervals, and Z = (Zi ,..., Z p )' is a vector of functions. The renewal theorem 
will give the limit of Z at +oo. 

For all real t, we set: 






TAIL OF MARKOV-SWITCHING DIFFUSIONS 


9 


(a) B = ( bij)i<i,j<p where b t j = J uFij(du) if it exists, the expectation of 


F 


(b) F^°\t) = where S t j(t) = l t >o if i = j and 0 otherwise, 

so that F-) * H = H for all H as in the definition above, 

(c) F^ n \t) = F * F^^it), the n-fold convolution of F, 

(d) U(t) = the renewal function associated with F. 

We will also assume that all the measures Fij are finite: 


Fij( oo) = lim Fij(t) < oo, 

t —KX3 


and that F( oo) is an irreducible matrix. F( oo) being an irreducible nonneg¬ 
ative matrix, we can use the Perron-Frobenius theorem: its spectral radius 
p(F( oo)) is a simple eigenvalue with right and left positive eigenvectors. We 
will also assume that p(F(oo)) = 1, and we choose two positive eigenvectors 
m and u so that 


p p 


F(oo)m = m, u , F(oo) = u', ^mj = 1, y] Ujmj = 1. 


We also assume that the sequence (||F(oo) n ||) is bounded [e.g., if F{ oo) 
is aperiodic, this is true]. We recall the following definition: F is lattice if 
the following conditions are satisfied: 

(a) For all i j, F,j is concentrated on a set of the form b t j + 

(b) For all i, Fa is concentrated on a set of the form AjjZ. 

(c) Each A a is an integral multiple of the same number. We take A to be 
the largest such number. 

(d) For all a^, a^, points of increase of Fjj , F t j ., respectively, 

a ij + djk ~ O'ik is an integral multiple of A. 

We can now state the renewal theorem. It extends a previous result of 
Crump (1970) and Athreya and Rama Murthy (1976) which deals with the 
positive case: each distribution F^ has support on R + . The proof of this 
theorem is given in de Saporta (2003). 

Renewal theorem A. Assume that F is as above and that, in addi¬ 
tion, it is a nonlattice matrix, that its expectation B exists, and that for all 
Rl, U(t) is finite. If G is directly Riemann integrable [see Feller (1966)], 
and Z = U * G exists, then for all i, we have 



where m and u are the eigenvectors defined above and c= ( u’Bm ) 1 (under 
these assumptions, v!Bm^D). 
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We also recall Theorem 2.3 of Athreya and Rama Murthy (1976) that will 
be used in Section 7. Note that this theorem can now be seen as a corollary 
of Theorem A. 


Renewal theorem B. Let F be a nonlattice matrix of distributions 
with support on the positive half-line, such that: 


(i) p(F( 0))<1, 

(ii) F( oo) is finite, irreducible and aperiodic, 

(iii) there exist i and j such that Fij{ 0) < Fij{ oo). 

Assume also that there is an a > 0 such that p(F a ) = 1, where (F a )ij = 
fo° e~ au Fij(du). Then for all h> 0, and all i,j, we have 


lim 

t—>oo 


rt-\-h 


e~ ay U i j(dy) 


cmiUjh, 


where rn and u are right and left eigenvectors of F a , with the same normal¬ 
ization as above, c= ( u'Bm)~ l , and B = bij with bij = f^°ue~ au Fij{du), c 
being interpreted as zero if some bij is equal to infinity. 


5. The renewal equations. Now we are going to derive the renewal equa¬ 
tions associated to our problem. In the following, we will suppose that the 
assumptions of Theorem 2 are satisfied. We set 5 = 1, and k will denote the 
unique positive solution of p(A s ) = 1. We are going to study the discretiza¬ 
tion yw. 


5.1. Notation. As A is a stationary process, we can extend it to negative 
t and define the coefficients <h n , V n and for negative values of n. Let 
bn = Vn 2 in and 

OO 

Bn — ) ^n^n —1 ' ' ' ^ n—k+l^n—k 

k =0 

(instead of Y n ) be the unique stationary solution of (4): R n+ 1 = 4> n+ i R n + 
b n . )_i. The limit law v of Y is also the law of R\. Thus we are going to study 
the random variable R\. 

The tail of the stationary solution of such recursive equations has already 
been studied in various cases. In the i.i.d. multidimensional case: <3? n are 
matrices and R n and b n vectors, renewal theory is used in Kesten (1973) 
to prove a heavy tail property when the <h n either have a density or are 
nonnegative. These results were extended in Le Page (1983) to a wider class 
of i.i.d. random matrices. Finally, in Goldie (1991) a new specific implicit 
renewal theorem is proved and the same results are derived in the i.i.d. 
one-dimensional case. This theorem also applies to the study of the tail of 
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several other random recurrences implying i.i.d. random variables. Recently, 
Goldie’s results were extended in de Saporta (2004) to the case where (<h n ) 
is a finite state space Markov chain. Here, (4> n ) is not a Markov chain, but 
conditionally to X n , <L n and 4>„ + i are independent. Our proof is thus very 
similar to that of de Saporta (2004), but we will repeat all the details for 
completeness. 

Note that £ n are standard Gaussian random variables, thus they are sym¬ 
metric, and they are also independent from the sequences (4> n ) and (V n ). 
Hence we have 


PJ >t) 

\k =0 / 

=p, J £ • • • ^2- k v!%- k > t) 

\k =0 / 

=p A E • • • ^2-kVlL 2 k (-^-k) > t 


v k=0 
' oo 


- ^2 • • ' d> 2-kbl-k > t 

\ k =0 ) 


Thus we have z/(]t, +oo[) = — oo,— 1[) for all t; hence if one of the limits 
stated in Theorem 2 exists, the other exists too and equals the same value. 
Therefore we need study only one limit. 

To study the tail of R \, we introduce a new function. For all t in M, we 
set 


z{t) = e 


-t 


u K ¥(Ri > u ) du. 


Lemma 9.3 of Goldie (1991) ascertains that if z(t) has a limit when t tends 
to infinity, then t K P(Ri > t ) also has the same limit. 

For all i in E and t in R, we also set 



u K P(R 1 >u,X 1 


i ) du, 


so that z(t ) = Y^,iLi Zi{t). We are now going to prove that Z = \Z \,..., Zjy) 
satisfies a system of renewal equations. 


5.2. The renewal equations. As R n satisfies (4), we have R\ = TiRo + ^i; 
thus for all t in R, we have 

P /i (Ri > u,X i = i ) = P At (d>i R 0 >u,X i = i) + 
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where 


ipi(t) =P M (i - b x < 4>ii? 0 < t, Ad = *) -P/i(i < 1 R() <t-bi,X 1 = i). 

We set Gi(t ) = e“ 4 / 0 e u K i/ji(u) du, and G = *(Gi,..., Gat). Then we have 

N r /*e £ 

z(f) = V e _< / rt K P„(d>ii?o > Xi = i) du + Gj(i) . 
i=i L J 

We have $1 > 0 and conditionally to Xq, $1 and i?o are independent. Thus, 
a simple change of variable and stationarity yield 


-t 


/■ e 

/ « ft P„( < I > ii?o > u, Ad = i) du 

Jo 

N ne t 

= ^e~ l / u K F fl ($iR 0 >u,X 1 =i\X 0 =j)n(j)du 

3=1 Jo 

N /»e £ 

= I u K ^j(^iRo >u,Xi = i)n(j)du 

T~, Jo 


3 = 1 

A r ^gt-log®! 

= E E i 

i=i L 

N r ggt-log®! 

= E E i 

3 =1 L 

iV r gt—logS^ 


tt K Pj(i?o > w) du 


»(j) 


tt K P M (i ?0 > w| A 0 = j) du 


A* O') 


= E E i *fl Xl =ie-( t - ] °s* 1 '> f u^^yu.X^fidu 

3 =1 L 7o 

Thus we get the following system of equations: for all i in E, we have 

N 

Zi(t) = J2 pEit^fl x 1= iZj(t - log^i)]] + Gi(t) 

3 =1 


(9) 


N 


= y £[Fij*Z j (t)) + G i (t), 
3 =1 


where Fij(t ) = Ej[$f lxi=ilt>iog$i]- Thus F = (Fij)ij^E is a matrix of dis¬ 
tributions in the sense of Section 4, and system (9) is a system of renewal 
equations that can be rewritten as Z = F * Z + G. To apply Theorem A, we 
now have to prove that F and G satisfy its assumptions. 
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6. Proof of Theorem 2, part I. As E is a finite set, <hi is bounded. There¬ 
fore, for all i,j in E , the measures Fij are finite and Ijj(oo) = Ej[3>ylx 1 = i]- 
Note that F(oo) = A' K . As A K is irreducible and aperiodic by Proposi¬ 
tion 2, so is F( oc), and its spectral radius also equals to 1. Besides, we 
have b t j = Ej[ < f>flx 1 =jlog^ ) i], thus the F\j have finite expectation. 

We are going to prove that the other assumptions of Theorem A are valid 
here in the following sections. 

6.1. F is nonlattice. Set a m = min.; G £{a(z)}, a« = nraxj e E{a(i)} and io, 
jo in E such that a(io) = a m and a (jo) = om- 

Proposition 7. For all i,j in E, x G }a m ,aM[ and small enough e > 0, 
we have 

Pj ^ a(X u ) du € }x - £]x + e[, Aj = jj >0, 
that is, x is a point of increase o/log<hi conditionally to Xo = i and Aj = j. 


Proof. Set x G ]a m ,aM[ and 0 < t < 1 such that x = ta m + (1 — t)aM- 
Fix i and j in E. As q is an irreducible matrix, we can find integers 0 < l < 
m < n and h,... , k n in E such that q^qkiM ''' <lk h i 0 > 0, qi 0 ,k l+1 qk l+1 k l+2 ■ ■ ■ qk, 
0 and qj 0 ,k m+1 qk m+1 k m+ 2 ■' ■ Qk n ,j > 0. Set also y = a{i) + a(h) + • • • + a(fcj) - 
(l + l)a m + 0(^+1) + ■ • ■ + a{k m ) — (n — l +1 )om +a(&m+i) + ■ ■■ + a(fc n ) + a(j), 
and z = min{£|y| _1 ,t(Z + l) -1 , (1 — t)(n — l + l)^ 1 }. Then we have 

P i(^J^ a(X u )due]x - e;x + e[, Xi=j S J 


>P i(X u = i on [ 0 ; 77[, X u = ki on [77; 2 r/[,..., X u = ki on [lip, (l + l)r/[, 
X u = i 0 on [(l + 1)77, t[, X u = k i+ 1 on [t; t + rj[, X u = k t+2 on 
( 10 ) [t + 77; t + 2 rj[,X u = km on [t + (m — l — 1)77; t + (m- l)rj[, 

X u = jo on [t + (m- l)rj ; l-(n-m + 1 )rj[, X u = k m+ ± on 
[1 — (n — m + l)?y; 1 — (n — m)r][,X u = k n on [1 - 2 77; 1 - rj[, 

X u = j on [1 — 77 ; 1]; r/e]0;z[). 

Indeed, on this event we have 

[ a(X u ) du 
Jo 

= rja(i) + rja(ki) H-h r]a(ki) + (t-(l + l)r})a m + rja{ki + i) 
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H-b ya(k m ) + ((1 -t)-(n-l + 1 )rj)a M + r]a{k m+ 1 ) 

H-b r]a(k n ) + ija(j) 


= ta m + (l-t)a M + VV 
= x + r]y, 


thus if y < e\y\ ^ then we have /g 1 a{X u ) du £ ]x — e; x + e[. Probability (10) 
can be computed [see, e.g., Norris (1998)]: 

(10) = fJ-(i)qi t kiQki,k2 ■ ' ' Qki,ioQio,ki +1 ' ' ’ Qkm,joQjo,km+l ' ' ' Qk n ,j 


x X(i)X(ki) ■ ■ ■ X{k n )X(io)(l - l)A(jo)(n -1 + 1) 


x 


/ [ e —A(*)r? e —A(fci)r? 

0 


• • • e 


A(fc n )?? e -A(io)(t-(/-l)» 7 ) e -A(jo)(l-t-(n-Z+l)» 7 e -A(j)r?j ^ 


Thus our choice of k \,..., k n and z ascertains that this probability is positive, 
which proves the proposition. □ 

Therefore none of the T)j(-) = Ej^Txj^T.Mog#!] can be concentrated 
on a lattice set, and in particular F is nonlattice. 

6.2. Finiteness of U. We are going to prove that for all i, j in E and t 
in K, Uij(t ) is finite. We start with computing the n-fold convolution of F. 

Lemma 1. For all n,i,j,t we have 



Proof. For n = 1, it is the definition of F. Suppose the formula is true 
for a fixed n. Then the Markov property and stationarity yield 




ft=i 



«=i 
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X [$£ +1 5 u (log $ n+ i) lx n+1 =i 1 X„=k] 


= Y 1 ■ • ■ ®Zhog$ 1 ---$ n <t-\og<l> n+ ilXo=j\lx n =k] 

ln(k) 


N 


k =1 


X E M [$£ +1 lx n+1 =i|lx n =fc] 




A 


= E m [$i ■ • ■ ^*n+lllog■ ■ -$nlog$„+! lx 0 =j 1 A„+i=i1 1 X n =fc 


k =1 


M(fc) 

™+l— 2 I A n = /cJ 


= E m [$i • • • ^^ + ili og $ 1 ...*„$„ +1 <tlx 0 =jlx n+ is 


1 


p(j) 


= Ej[$f • • • ^^ +1 liog*i-®„$„+i<tlx„ +1 =i]- 
Thus the formula is also true for n + 1 and the lemma is proved. □ 

We have seen that F( oo) = A Proposition 1 and the preceding lemma 
also imply that for all n we have F^ n \ oo) = (A™)' = F(oo) n . We can prove 
a more general result. 

Lemma 2. For all n and 0 < r < k we have 

/ OO 

e- ru F^ n \du) = {A n K _ r )'. 

-OO 

Proof. For all i , j in E, Proposition 1 and the preceding lemma yield 
e-W^Mu) 

D 

/ OO 

e-"%-[$f • • • $£<5„(log • • • $n)l^=«] 

-oo 

= E i [$f---^e- rlo «* 1 -*"l Xn =i] 

=E J [d»r r ---$r r iA„= i ] 

= (4U)ii- a 


Now fixO<r<ft. We have 

oo 

Ml 


u^t) = Yw^< 


„rt 


n =0 


t 00 

-r« \ 1 z?( n ) / 




n=0 


< e r ‘ E / e ~" F if (*•) =a E «-A<. 


n=0 “ 


n=0 


(ii) 
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according to the preceding lemma. But Corollary 3 says that p{A Kj ^ r ) < 1. 
Thus the series in (11) converges. Hence Uijft ) < oo for all i . j in E and t in 


K 


6.3. Proof of Z = U * G. Iterating the renewal equation (9) yields, for 
all n, 


n —1 


Z = F (n) * z + V F (fc) * G. 


( 12 ) 


k =0 


The same change of variable as in Section 5.2 yields 



i= 1 


But we have seen at (7) that we have ‘hi • • • —> 0 when n tends to infinity. 
Thus the bounded convergence Theorem yields J2iLi{F^ * Z)i(x,t ) —> 0 as 
n tends to infinity. Each term of this sum is nonnegative, thus each term 
tends to 0. Letting n tend to infinity in (12), we thus get Z = U *G. 

6.4. G is directly Riemann integrable. As the Gi are clearly continuous 
in t, it is sufficient to prove that 


OO 


Y sup |Gj(f)|<oo 

l=-oo l < t<l + 1 


[see Feller (1966)]. But for all i,t , we have Gift) = G}(t) — Gj(t), where 




For all real t, we have Gift) < G}(t) < e f fff u K du = e *( re + 1 )( K -g i) 1 , i n 
particular, Gi is directly Riemann integrable on K_. We still have to study 


G\ and Grf on M + . These two functions being of the same kind, we only give 


the detailed study of the first one. 

The proof is adapted from Le Page (1983). Set e € ]0; 1[ such that —1 < 
k — (1 — e) < 0. Thus we have 



(13) 



+ / u h P^(u — u £ < ‘hiRo < u, X\ = i) du. 


We are going to give an upper bound for each one of these two terms. 
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First term. Chebychev inequality yields 



1 + /c(l — e) ’ 


,t(l+s(l-e)) 


Note that b\ has moments of all order. Indeed, we have, by independence, 
!E^|&i| K = IE//.(I / 1 K ' /2 )IE; t |^i| re , and £1 is a standard Gaussian variable and V\ is 
bounded. 

Second term. We have 





< [ u K [\ — l u >i(u — u £ ) K (l — eu £ 1 )\§ D n(QiRo>u — u e ,X$ = i)du. 
Jo 

Set 0 < r < At. As <f>i is bounded, there exists a positive constant c such that 
for all u > 0 we have 



which is bounded by Proposition 6. Thus we get 

re 1 

(15) / u K FJu-u £ <^> 1 R 0 <u,X 1 =i)du<Ce t ^- r+e - 1 \ 

Jo 

where C is a positive constant. Now set (3 = ma x{k + £ —r; l + /c — ks} € ]0; 1[. 
Then (13)-(15) yield e t G\ (t) < ce t J for all t > 0. Thus Gj(t) < is 

directly Riemann integrable on M + . 

6.5. Tail of the distribution. We have now proved that F and G satisfy 
the assumptions of Theorem A. Thus we get, for all i,t , 



N 


•oo 


(16) 


Gj(y)dy. 


— OO 


Summing up these terms, we get 


(17) 



as Yl m i = I- We still have to prove that this limit is nonzero. 
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7. Proof of Theorem 2, part II. Now we are going to prove that there 
exists a positive constant C such that t K P M (|Ri | > t) > C > 0 when t tends to 
infinity. First, we give a lower bound of this probability involving the prod¬ 
ucts 4>i • • • 4> n , and then we study the asymptotic behavior of such products. 

7.1. Lower bound for v{x : |x| > t}. The following proof is adapted from 
Goldie (1991) and de Saporta (2004). 

Proposition 8 . There exist e > 0 and a corresponding positive constant 
C such that for large enough t we have 

M^ll^^^SUp ($!•■■$„) >|). 

For the i.i.d. case, the key to such a lower bound is an inequality estab¬ 
lished in Grincevicius (1980) that extends Levy’s symmetrization inequality 
[see Chow and Teicher (1978)]. Here we need first to extend this inequality. 

Recall that R\ = Y^kLo • ■ ■ ^ 2 -kW-k- For all n > 1, we set 

n— 1 

R\ = $>i<&o • • • 3*2—fcfoi —k an d n n = <I>i<I>o • • • ®2-n- 

k =0 

If x is a a(X tl Wt,a < t < immeasurable random variable, let medj(a:) be a 
median of x conditionally to = i and med_(x) = minj{medj(x)}. 

Lemma 3. For all t > 0 and n > 1, we have 

P M ^ max + n, med_ | > tj < 2P M (R” > t). 

Proof. Set T = inf { j < n t.q. R{ + med_(nF 1 (R” — Rjj) > tj if 
this set is not empty, n + 1 otherwise, and Bj = {med_(nj 1 (i?” — Rj)) < 
n J 1 (Rf — R \)}. The event (T = j) is in the cr-field generated by (X t ,W t , (1 — 
j) < t < 1), and Bj is in the cr-field generated by (Xt,Wt, (1 — n) <t < 
(1 — j)). Therefore these events are independent conditionally to Xy_jy Be¬ 
sides, for all i and j we have ¥ fl (Bj\X^ 1 _j' ) = i) > P /i (medj(nF 1 (R^ — R\)) < 

nj 1 (R™ — R\)\X(i_j) = i) > 1/2. Thus, as the products Hj are positive or 
zero, we have 

P M (R?>f)>Pj (j[B j n(T = j)} 

\j=i 
n N 

= E E = iMT = j |X (1 _ j} = i)p(i) 

3 =1 *=1 
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> )fn(T < n) 

=MH + u> med i^r) } > *) • □ 

Under our assumptions, R™ tends to R± when n tends to infinity, and 
for fixed j, — R{) converges to a random variable R that has the 

same distribution as R±. Set mo = med_(i?i) = min,;{med(/?,i \X\ = i)} = 
med_(-R), and letting n tend to infinity in Lemma 3, we get, for all t > 0, 

^sup{/?,"[ + IL,m 0 } > tj < 2F fl (Ri >t). 

Replacing R\ by —R\ yields a similar formula; thus, for all t > 0 we get 

(18) ^V^sup \R{ +U j m 0 \ >tj < 2P m (|Ri| > t). 

Furthermore, as proved in Goldie [(1991), page 157], for all t > |mo| we have 

P H ^ sup {R™ + U n m 0 } > tj 

> P M (3n s.t. \(Ri +1 + n n+ im 0 ) - (i?i + n n m 0 )| > 27), 
where R® = 0 and IIo = 1. But we have 

(Ri +1 + n n+1 m 0 ) - (7?i + n n m 0 ) 

= •Fld’o ' ' ' ^2-nbl-n + (Lln+1 — Ll n )mo 

= n n (f>i_ n + (d»i_ n - l)m 0 ). 

Thus (18) yields, for all e > 0, 

P^d-Ril > t) > ^P M (3n s.t. |n n (6i_ n + (<Fi_„ - l)m 0 )| > 2 1) 

(19) 

> 

Now we give an extension of Feller-Chung’s inequality adapted to the 
present case [see Chow and Teicher (1978)]: 


ip^EIn s.t. |n n | > j and |6i - n + (3>i_ n - l)m 0 | > sj . 


Lemma 4. For all t > |mo| and s > 0, we have 


P s . t . |n n | > j and \bi- n + (®±- n - l)m 0 \ > sj 
> min Pj(|6 0 + (4> 0 - l)m 0 | > e)P s . t . \ U n \ 

l<.i<.N \ 
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Proof. Set A 0 = 0, A n = {|n n | > 2te~ 1 } and B n = {\b\- n + (4>i_ n - 
l)mo| > s}- Conditionally to Xn_ n ), B n is independent of Aq, ... ,A n . Thus 
we have 


IP J U [ A n n B n 


k n=l 


oo 


71—1 


= ^2 ( B n n A n fi [Bj n Aj] 

j =0 


n =1 
oo 


71—1 


> £ (B n n An PI A c j 


71—1 
oo N 

= EE 

n =1 i =1 


3=0 


71—1 


^li(B n \X(i _ n ) — i)P^ f'j n) — m(®) 


where H c denotes the complementary set of A. But, by stationarity we have 
P M (S n |X (1 _ n) = i) = Pj(|6 0 + (<f> 0 - l)m 0 | > £)■ Thus we get 

OO \ / OO 

(J [A n n B n \ J > min Pj(|6 0 + (<f> 0 - l)m 0 | > e)P M ( |J A n 
n= 1 J X < l < N \n=l 




Now we can give the proof of Proposition 8. 


Proof of Proposition 8. Equation (19) and Lemma 4 yield, for all 
t > |mo| and e > 0, 

P M (|i?i| > t) > - min Pi(|6 0 + (d> 0 - l)m 0 | > e)P M fdn s.t. |II n _i|>—Y 
Z l<i<lv \ £ J 

1/2 

We have bo = V 0 £o, Vo and 4*0 are bounded, but £ is not bounded as it is a 
Gaussian variable. Thus equality b o + (4>o — l)mo = 0 cannot hold P,;-almost 
surely. Thus we can find e > 0 such that mirii<j<jv P,;(|6o + (4*0 — l)mo| > 
e) > 0. Hence, as expected there is a constant C > 0 such that for all t > |mo|, 
we have 

P/i(|fli| >t) > CP M ^sup |n n | > n 

7.2. Asymptotic behavior of the products dq • • • 4> n . To estimate the prob¬ 
ability P M (sup„ |n n | >t), we use the ladder height method given by Feller 
(1966) for the study of the maximum of random walks. 
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7.2.1. Notation. First we introduce some notation. Set So 
all positive n, we set 


s n = V log(<F 2 -fc) = logn n = / a(X u ) du. 

J C 1 -") 


0 and for 


The first ladder epoch of this random walk is r = n = inf{n > 1 s.t. > 
0}, and the first ladder height is S T . We denote by H(t) the matrix of 
distributions of S T with the following coordinates: 


Hij(t) = P m (t <oo,S t < t,X( i_ T ) =j\Xi =i). 

As S T > 0, H is distributed on the positive half-line. Moreover, S T > 0, 
Si-r < 0 and the <F n are bounded, thus we have SV < sup log < sup a{i) < 
oo, and H has bounded support. 

We define also the nth ladder epoch by r n = inf{n > r n _i s.t. S n > S Tn _ 1 }, 
and S Tn is the corresponding ladder height. We check that we have 

( t) = P M (r n < oo, S Tn < t, X {l _ Tn) = j\Xi = i), 

where H^ is the n-fold convolution of H. Let 'F = be the renewal 

function associated with H. 


7.2.2. The random walk S Tri . To investigate the asymptotic behavior of 
(S Tn ) we are going to use a renewal theorem as in Feller (1966) for the i.i.d. 
case, namely, Theorem B. We want to apply it for F = H and a = n, thus 
we have to prove that H satisfies its assumptions. 

As H( 0) =0, we have p{H{ 0)) < 1, thus the first assumption is true. In 
addition, H t] are probability measures, therefore H is finite. H has bounded 
support because S T - 1 < 0, S T > 0 and $ is bounded. Thus B, the expectation 
of H k { oo) = / 0 °° e~ KU H(du) is well defined. Proposition 7 yields again that 
H is also nonlattice. 


Irreducibility and aperiodicity. For all i.j in E, we have 
Hij(oo) =P M (r < oo, X 1 _ T =j\X l = i) 

>P M (r = l, X 0 =j\X 1 = i) 

— Pj(log d>i > 0, X 1= i)^4 

p{i) 

= P J ( Jo a (Xu) du > 0, X, 

and Proposition 7 implies that the last term is positive as 0 G ]a m ; om[. Thus 
the second assumption of Theorem B is valid. We have also proved that for 
all i and j we have 0 = Hij( 0) < Hij(oo), so that the third assumption is 
also valid. 
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Spectral radius of H K (oo). Now we define a new probability law P K on 
11x0. For all bounded A x ^-measurable functions / which first coordinate 
depends only on (X t , (1 — n) < t < 1), we set 


P«(/)=E«(/) 

_ E M (/(ch 1 ,...,«l> 2 _ n ,g)(cl) 1 ...«h 2 _ n ) K ) 

E M (($i---$ 2 _ n )«) 

Set H K (t) = Jq e~ KU H(du). We have 


(HMt) 


E k (t < oo, S r < t, X (1 _ r) = j|Xi = *) 
E At (($i • • •d>i_ r ) K ,r < oo) 


E/i((d>i • • • d>i_ T ) K , r < oo) ’ 

where = P K (r < oo, S T < t, ^(i- T ) = j\X\ = z) describes the behav¬ 

ior of the ladder heights of our random walk under the new probability 
law P re . 

The computation we made in the proof of Proposition 3 yields 

log*h* 


d_ 

dr 


log 00^))= lim —Ej 

n—>oc 77, ' zJ 


\ i=l 


= E K (log$i). 


But we have logp^o) = logp(^4 re ) = 0; this function is convex (Corollary 
2) and its right-hand derivative at 0 is negative (Proposition 3). Thus its 
left-hand derivative at n is positive, that is, E h ,(log<l>i) > 0. Under the law 
P K our random walk thus drifts to +oo, hence for all n and i, we have 
(P K ),(r n < oo) = 1 and is a stochastic matrix, therefore its spectral radius 
equals to 1. 

For all n, we have 


tf<”>(oo) = (ff„(oo))" 


(gKoor 

E At ((<f>i • • • <f> 2 _ T J K , r < oo) ’ 


thus p(H K ( oo)) = lim(E At ((d>i • • • <h 2 _ T7i ) K , t < oo)) -1 /” and we now have to 
prove that this limit equals to 1. But for all n, we have r n > n, and the event 
(r n = k ) depends only on ( X t , (1 — k) < t < 1). Thus we have 


E M ((<Jh T n <00) 


= YM^i--^i-k) K ,r n = k) 

k=n 


= ^P K (r n = fe)E ft (($ 1 ---$ 1 _ fc ) K ). 

k=n 


(20) 







TAIL OF MARKOV-SWITCHING DIFFUSIONS 


23 


Set e > 0. For large enough n, our choice of k and (5) and (6) yield 
M" 1 - e < ■ ■ ■ t’l-nT) < nA n K 1 + e. 

Thus for large enough to, (20) yields 

OO 

(M"l - <0 E P «( T « = k) 

k=n 

oo 

< E m (($i • • -$l - Tn ) K ,T n < Oo) < 1 + e) F K (r n = k), 

k=n 

and as P K (r„ < oo) = 1, we have 

M” 1 - £ ' ' ' ^l-r n ) K , T n < OO) < IXA™ 1 + £. 

Thus as n —> oo we have, with the notation of Corollary 1, E^(3>i • ■ •<]? 

1. Hence we have, as expected, E m ((3>i • • • <bi_ Tri ) K ,r n < oo) 1//n —> 1. 
Thus all the assumptions of Theorem B are valid here. We are going to 
use it in the following part. 


7.2.3. Asymptotic behavior of the maximum. Let M = sup n S n = sup n S Tn 
be the maximum of our random walk. Using the definition of H, we get, for 
all 1 < i < N : 

F li (M<t\X 1 =i) 

OO 

= Y F //( r n < °°> S Tn < L T n+ 1 = Oo|Xl = i) 
n= 1 

oo N / -\ 

< OO ,S Tn < t,T n+ 1 = oo,Xi = z|X (1 _ rn ) =j)~J- 

n =1 j=l 


oo AT 

( 21 ) = X] E [ P f( t « < °°> 5 Tn < t, X {1 _ Tn) = j\X, = i) 

n=l j=l 

x (1 - P /i (r n+ 1 < oo\X (1 _ Tri) = j)} 

oo N 

= EE 

n=l i=l 
A 

=E 

3= 1 


( 1 - E 


A 


fc=l 


A 


¥)i-E ff i‘( 

V fe=i 


oo 


Theorem B applied to (21) yields, when f tends to infinity, 
l-P /i (M<t|X 1 =z) 
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( 22 ) 


N 

=E 

j =1 L V fc=l 


t—> oo 

rvj 


1-E^Mj r^e^du) 

k=\ ) Jt 

‘ / \ /*oo 

E — y~) iJ ? fc(oo)J J e~ KU c7niUjdu 


N 


= E 


AT 


1 - E H jk (oo) ) cfhiUj 


j=1 L \ fc=l 


— Kt 


where m and u are right and left eigenvectors of H K (o o) with positive coor¬ 
dinates with the same normalization as in Section 4, and c= (* uBfh ) _1 > 0. 


7.3. Conclusion. We still have to prove that there is a j < A r such that 
1 — J2k=i Hj k (oo) > 0. But the mapping r i—+ H r (oo) = f^°e ru H(du) is 
clearly increasing component-wise. As these matrices are nonnegative and 
irreducible, Corollaries 8.1.19 and 8.1.20 of Horn and Johnson (1985) im¬ 
ply that the mapping r i—> p(H r (oo)) is also increasing. As p(H K (oo)) = 1, 
we have p(Hq(oo)) = p(H( oo)) < 1. This is a substochastic, nonstochastic 
matrix, thus there exists a j such that we have 1 — J2k=i ^jfc(°°) > 0. 

We have now proved that the right-hand side term in (22) is positive, thus 
there is a constant C > 0 such that, when t tends to infinity, we have 

N 

(23) e K ¥ M (M >t)>J2 e Kt r^(M > t\X x = i)p{i) > C. 

i— 1 

Putting together this result and Proposition 8, we get, for large enough t, 

(24) t K P M (|i?!| >t) > K>0. 

With the notation of Theorem 2, it means that L >0, which ends the proof 
of this theorem. 


8. Determination of k. Set si = min{A(i)a(i) _1 |a(i) > 0}, and let M s 
be the matrix with components {q(i,j)\(i)(\(i) — sa(z)) -1 }. This matrix is 
well defined for all s < si. We can precisely compare the spectral radius of 
A s and that of M s . 

Proposition 9. For all 0 < s < s\, we have p{M s ) < 1 if and only if 
p(A s ) < 1, and we have p(M s ) > 1 if and only if p(A s ) > 1. 

Proof. Suppose that p{M s ) < 1. M s is a positive irreducible matrix 
as q is, A being positive and s < si. Thus the Perron-Frobenius theorem 
[see, e.g., Horn and Johnson (1985)] gives the existence of a vector tp with 
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positive coordinates such that M s ip = p(M s )tp < p. Hence for all i in E, we 
have 

<p{i) 

j A (z)-sa(z) 

that we can rewrite, since s < si, as 

(25) (sa(i) - A(i))<p(i) -h < 0 . 

3 

Proposition 4 enables us to choose a small enough 5 such that (8) is valid 
here. Equation (25) thus yields 

A s p(i) = [1 + 5{sa(i) - A(i))]<p(i) + 6X(i) ^[<?(i, j)<p(j)] + o(S) 


= tp(i) + 5 


(sa(i) - X(i))ip(i) + \(i)^q{i,j)y(j) 

3 


+ o(5) 


< ¥>(*)■ 


Thus component-wise we get A s p < p, which implies that p(A s ) < 1. The 
proof that p(M s ) > 1 implies p{A s ) > 1 runs the same, the inequalities being 
reversed. 

Suppose now that p(A s ) < 1. A s is a positive irreducible matrix, thus the 
Perron-Frobenius theorem gives the existence of a vector i/j with positive 
coordinates such that A s if = p(A s )ip < i/j. Hence for all i in E, and small 
enough 5, we have 


5 


(sa(i) - A (i))i/>(i) + A(i) J2 ?(*> j)V’(j) 

3 


+ o(6 ) = A s ^(i) - ipi 


< 0 . 


Hence, for all i, we get ( sa(i ) — X(i))ijj(i) + X(i) Yhj q(i, j)ip{j) < 0, or > as 
s < si, 

^(0 > 

A(«) - sa(i) j 

and thus M s if < if. As M s is a positive matrix, we conclude that p(M s ) < 1. 
Here again the proof that p(A s ) > 1 implies p(M s ) > 1 runs the same with 
reversed inequalities. □ 


Proposition 10. The spectral radius of M s tends to infinity when s 
tends to sq. 








26 


B. DE SAPORTA AND J.-F. YAO 


Proof. Set io £ E such that A(*o)a(*o) 1 = si, and e^ the row vec¬ 
tor with zero coordinates except the *oth which is set to be 1. Set Vi 0 = 
A(z 0 )(A(i 0 ) - sa(i 0 )) _1 . We have e io M s = v io q(i 0 ,-) > v io e io as q is a pos¬ 
itive matrix. As M s is also positive, for all s < si, we get p{M s ) > Vi 0 = 
A(zo)(A(zo) — sa(io))^ 1 - Hence this spectral radius tends to infinity when s 
tends to si. □ 


Corollary 4. There is a unique s £ ]0;si[ such that p(M s ) = 1, and 
this s equals the unique k such that p(A K ) = 1. 

Proof. For all s < k, we have p(A s ) < 1 by Corollary 3; thus Propo¬ 
sition 9 yields p{M s ) < 1 for all 0 < s < min{«;, si}. As p(M s ) —>• oo as s 
tends to si, we also have p{A s ) > 1 for s close to ,s'i. Therefore k < s i, and 
p(A s ) > 1 for all k < s < si. Hence p(M s ) > 1 for all k < s < si. As M s has 
continuous coordinates, its spectral radius is also continuous; thus p(M K ) = 1 
and k is the only value of s £ ]0;si[ satisfying this equation. □ 


We now give an illustration by computing the value of k when E = {1,2}. 
The jump kernel q then equals to 


Q = 


0 1 
1 0 


and the invariant law of the process X is p = (A(2), A(1))/(A(1) + A(2)). We 
suppose that a(l) or a( 2) is positive. Condition 2 becomes 


(26) A(l)a(2) + A(2)a(l) < 0. 

For all i in E, set r* = We have rq + r 2 < 0, r\r 2 > 0 and si = rnaxjrf 1 , r^ 1 
For s £ [0, si [, the matrix M s equals to 


M s = x 

V 1 - sr 2 

and its spectral radius is [(1 — S7q)(l — sr 2 )] _1 - /2 . It equals to 1 for k = 
rf 1 + r^ 1 = A(2)a(2) _1 + A(l)a(l) -1 . 
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